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A formal formulation of the differential cross section for x-ray inelastic scattering is given for a 
real solid, in particular, in terms of the polarization propagator and the inverse 1 dielectric function. 
The differential cross section is related to the causal functions of electron properties rather than those 
retarded functions. 
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A few years ago, the present author formulated the 
differential eross section for x-ray inelastic scattering 
in a unified way such that core electrons as well as 
valence electrons were treated on an equal basis [1J. 1 
This formulation aided explanation of the then un- 
resolved x-ray Raman scattering problem. In this 
paper, tbe differential cross section will be expressed 
in terms of causal rather than retarded functions in 
model-independent form. The results are similar in 
form to the Nozieres-Pines [2] results for an electron 
gas in the dielectric formulation. 

Several results of this paper have been reported 
previously or used separately in papers published 
elsewhere [3, 4]. It is, however, desirable to present 
the entire formulation in a systematic way, although 
the formulation is a trivial one to those who are familiar 
with Schwinger's many-body formulation and the spec- 
tral function technique. 

In previous papers [1, 5] it was shown that the dif- 
ferential cross section for x-ray inelastic scattering is 
given by the double four-dimensional Fourier transform 
of the electron charge (or current) correlation: 
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where k and k' are the four-vectors consisting of the 
wave vectors and the energies of the incoming and the 
outgoing photons, respectively, and € and e' are their 
polarization directions. The quantity C(k, k') is defined 
by 



C(k, k')=(2TT)- 4 {e 2 lm) 2 (kk')-K 
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The function p(x) is the electron charge operator and 



Ap(x) 



pven by 



Ap(*)=p(*)-<p(*)>, 
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where (Q) indicates the expectation value of an 
operator Q with respect to the ground (Fermi vacuum) 
state of the electron system. This ground state 2 is a 
Heisenberg state of the electron system in the absence 
of radiation fields flj. Units are such that h = c=\ 
and the metric is kx — k ■ r— |Z :\t if the four vectors 
are given by k= (k, \k\) and x— (r,t) . 

Since the charge correlation in eq (1) is defined 
by the Fermi vacuum expectation value of the system 
of electrons per se, one can proceed via the standard 
approach for many interacting electron problems to 
calculate it. In particular, we use Schwinger's action 
principle [6] by adding an external c-number source 
function to the Hamiltonian density of the electron 
system. The Hamiltonian density due to the external 
source, U (x) , is given by 



f s (x) = U(x)^(x)il>(x), 



(3) 



where i// (x) is the electron field operator and i// + (x) 
is its adjoint. Then the action principle gives 

i 8 8U(2) ) = ( 7, [<?(Dr(2)^(2)]> 



<<?(!)> <r(2)<M2)> 



(4) 



1 Figures in brackets indicate the literature references at the end of this paper. 

2 At finite temperatures, the ground state average (...) must be replaced by an averag 
over a grand canonical ensemble. 
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where the variable, say 1, is short for «i = (ri, £1); 
Q (x) is any operator in the Heisenberg representation; 
the symbol T [. . . .] implies the time ordering opera- 
tion for fermions. The electron Green's function is 
defined by 



G(12) = (-i) <r[*(D* + (2)]>. 



(5) 



From this, the vacuum expectation value of the electron 
charge operator, p (x ) , is given by 



<p(l)> = (fd + WD) =(-i)G(i,i + ), 

where 1+ indicates the four vector (ri, t\ + 0). 
If we substitute (6) into (4), we obtain 



(6) 



iS (12) = i 



§ <p(l)> = 8G(l,l+) = 
817(2) 8f/(2) 
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This quantity is a causal Green's function, while the 
quantity in eq (1) is merely a time correlation. These 
two are, of course, different functions. However, 
there is a definite mathematical relation between' 
them in the spectral representation [7, 8], provided 
that translational time invariance holds. In our electron 
system, this invariance holds since the total Hamilton- 
ian of the electron system is thought of as being time 
independent. 

At absolute zero, 3 the spectral relationship between 
the time correlation and the causal function is given by 

<Ap(l)Ap(2)> w = -2Im[(-i)<nAp(l)Ap(2)]> tt ], 

(8) 

where (. . .^ implies the spectral intensity (the 
Fourier transform using a basis of exp [ — ia> (ti — h) ] ) 
of the function (. . .), and the symbol Im denotes 
"the imaginary part of". Since the differential cross 
section, eq (1), contains the quantity given by eq 
(8), we now study the properties of the causal function 
eq (7). 

The electron system which we are dealing with must 
represent that of a real crystal. In such a system 
translational invariance may not hold for spatial 
coordinates. Therefore we deal with quantities, 
such as eq (7), as they are written, in terms of space- 
time coordinates; but not in terms of their Fourier 
transform until we reach the very end. The electron 
Green's function, G(12), satisfies the electron Green's 
function equation: 

G 1 (lT)G(l2)=G(lT)G- 1 (12)=8(12), (9) 

where the convention for integrations applies to the 
repeated variables with a bar. The vertex function for 
electron-electron interaction (not a radiation part) 
is defined by 



The quantity G -1 can be given in terms of the total 
average potential, F(l), of the system and the self- 
energy operator 2(12): 
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where R„ denotes the position of the nucleus n, 
Z n \e\ is the charge of the nth nucleus, v c is the instan- 
taneous Coulomb interaction, and Wj is the Fermi 
energy. The self energy operator X is given by 



2 (12) = -it; c (ll)G(12)rV(22; 1 + ). 



(13) 



The time ordered product of Ap has been connected 
to the variational derivative of the electron Green's 
function by eq (7). In terms of the vertex function, 
eq (10), we can write 



S(12) = iG(ll)Ml 2;2)C(21+). 



(14) 



This expression provides a way of calculating the 
charge correlation from eq (7). However, there is 
another way [9] to express it, emphasizing the dynamic 
properties of interacting electrons. Such dynamic 
properties of electrons can be characterized by a 
generalized dielectric or inverse dielectric function 
of the system, which represents the true response of 
the system, including the dynamic screening effect 
of electron interactions. The inverse dielectric function 
can be defined by 



(12) = 



8V(1) 
St/ (2) 



Then we obtain 



S(12)=P(ll)r'(12), 



(15) 



(16) 



where P is the irreducible polarization propagator 
[9, 10] . This propagator is now given in terms of a new 
dressed vertex function T by 



P(12) = iC(ll)r(12;2)G(21 + ; 



where 



r(12;3) = 



SG- 1 (12) 
8V(3) 



(17) 



(18) 



3 At finite temperatures, the right-hand side of eq (8) should he multiplied hy the following 
factor 



f(<oJ) = 



(8') 



I- 17 exp [ — oilkT] ' 
where T is absolute temperature, k is Boltzmann's constant and 77 is — 1 for fermions 
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We have thus calculated the causal function S(12) 
originally denned by eq (7). We have also found that 
the spectral intensity of the time correlation is related 
to that of the causal function by relation (8). Since the 
differential cross section (l _ a) is given by the Fourier 
transform of the electron charge time-correlation, it 
is now possible to calculate the differential cross 
section (1-a) from the causal function S(12) through 
eq (16). We write the double Fourier transform of 
function F (12) as follows: 



F(12) 



h 



do) dki c?k2 



2tt (2tt) 3 (2tt) : 



F (ki,kt;aO 



(19) 



exp [ikiri — ik 2 r 2 — iw (ti — t? ) ] . 



In terms of these Fourier transforms, eq (8) gives 
N(k u k 2 ; co)=-2 ImS(ki, k 2 ; <y), 



(20) 



where /V(ki, k 2 ; a> ) is the Fourier transform of the 
charge correlation (Ap(jc)Ap(y)). It is convenient 
at this stage to introduce two experimental variables. 
One is the momentum transfer k which is given by 
both the energy (or wave-length X in ) of the incoming 
x rays and the scattering angle 0. The other is the 
energy transfer E which is the energy loss of the out- 
going x rays. These quantities are defined by 



\k\ - (47r/\ in ) sin (612) 



and 



E 



(21) 



(22) 



Then the differential cross section (1-a) is rewritten 
in terms of the double Fourier transform of the charge 
correlation, N(k, k, E) as 



(r(k-*k')=(r(K 9 E) = 

C{k,k' )(€•€' 



'J* 



N I 



k, k: 



(23) 



It should be noted that the differential cross section 
defined by (1-a) is the total transition probability 
instead of the customary transition probability per 
unit time (i.e., rate). As seen from (23), the customary 
transition rate a (k, E) is given by the integrand of eq 
(23) since the integrand turns out to be time 
independent: r 

a(K,E) = J dta(K,E). (24) 

Finally the customary differential cross section is given 
via eq (20) by 



(T(k,E)=(-2) C(kk') (e-€') 2 



Im 



f$pPi*.*E) 



c- ! (k, k;E] 



(25) 



where the Fourier transform of eq (16) is used to ex- 
press S(k, #c, E). This expression has been used to 
calculate the inelastic scattering profiles at different 
scattering angles in both the Hartree-Fock and the 
random phase approximations [3, 4| . 

As mentioned previously, the present formulation is 
model-independent and also holds for real crystals. 
In such crystals, periodic spatial translations! invar- 
iance may not hold, thus forbidding the use of a Fourier 
series-integral expansion. As seen in the expression of 
the differential cross section, the final result is given 
by the double Fourier transforms of the physical quan- 
tities. Even in the crystal with imperfect spatial trans- 
lational invariance, one can always expand the nonlocal 
functions by the double Fourier integrals as we have 
shown here. 

Another advantage of the present result is found in 
its capability of separating, for instance, the core elec- 
tron contribution from the valence electron contribu- 
tion, thus filling in the gap between the modern x-ray 
inelastic theory and the traditional Waller-Hartree 
expression [11]. This separation can be done, because 
the polarization propagator, P, can be expressed as a 
sum of two terms: P = P CO re+ ^valence* It * s a ^ S(> possible 
to deal with the interaction of electrons in crystals 
realistically through a proper decomposition of the P: 
for instance, P = P d iagonal + Pnondlagonal- 

As a final point, note that, substituting expression 
(12) into (15), we obtain 



e" 1 (12) = 8 (12) + i; c (11)5(12). 



(15') 



The Coulomb potential, i; r (12), is instantaneous and 
a function of 1-2. Therefore, the Fourier transform 
of v c can be diagonalized. Utilizing this property, one 
may express the imaginary part of 5(12) by that of 
e _1 (12) , like in an electron gas model. 
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